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ABSTRACT 

The  mode  parameters  for  GPS  common  view  time  transfer  is  estimated  and  an  unbiased 
weighting  method  (UBWA)  for  data  in  GPS  common  view  is  discussed  in  this  paper.  The 
one  SI  jitter  of  each  time  series  is  determined  by  the  standard  deviation  of  residuals  of 
polynomialregression,  which  define  weights  used  to  compute  a  weighted  average.  Averaging 
the  double  difference  observables  during  an  integral  length  (3 0  days),  the  systematic  biases 
existing  between  every  two  time  series  are  corrected  in  constant  monthly  relative  to  a  reference 
time  series.  The  GPS  common  view  data  among  NRC  in  Ottawa,  NBS  in  Boulder  and 
USNO  in  Washington  DC,  are  calculated  by  using  three  different  methods:  UBWA,  equal 
weights  average  method  (EWA),  and  direct  weighting  method  (DWA).  The  final  weighted 
averages  are  compared  in  terms  of  T.  (standard  deviation  of  residual),  Allan  variance 


Results  showed  that  UBWA  is  superior  to  the  others. 


ESTIMATE  FOR  MODE  PARAMETERS 

GPS  single  difference  observables  between  stations  reflect  receiver  clock  differences  between  stations 
directly  and  apply  to  high  precision  time  comparison  widely,  i.e.  the  so  called  common  view  mode.  Its 
model  equation  can  be  expressed  as; 

l  /c{R  BJ(x))  +  «/ „-  f  j)/ f  r-9  f  A  -  f ,)/ f  fix  -  CQ)+  t  AB(ion)+  1  A[)(trop)  (1  ) 
where  A  and  B  refer  two  stations,  i  refers  the  satellite,  t  is  the  epoch  for  measurement.  lAJ(x)  an  d  t  BJ  ( i; ) 
are  the  time  readings  at  stations  A  and  B  respectively,  and  the  corrections  have  been  made  on  these 
readings  for  the  offset  of  the  receiver  system  time  delay.  R  ,,,/and  R  Bj  are  the  distances  from  satellite  j 
to  stations  A  and  B  respectively,  f  A  and  f  B  are  the  receiver  clock  frequencies  at  stations  A  and  B 
respectively,  f  t  is  the  satellite  clock  frequency,  t  AH(ion)  and  l.Ali(lrop)  arc  the  differenced  iono¬ 
sphere  error  correction  and  the  differenced  troposphere  error  correction,  respectively. 

A  common  view  track  which  repeats  every  sidereal  day  defines  a  time  series  comparing  the  clocks  at 
the  two  stations.  Each  time  scries  represents  the  time  measurement  of  signal  propagation  over  a  specific 
path  from  one  station  to  the  satellite  to  another  station.  Change  in  path  is  caused  by  different  satellite 
or  different  measurement  time,  and  forms  different  time  series.  Some  experiments  have  reported  that 
the  different  time  series  between  two  stations  has  not  only  the  different  time  jitter,  but  also  obvious 
systematic  bias.  Through  analyses  of  GPS  various  difference  observables,  the  reason  for  existing  jitters 
and  biases  could  be  understood.  For  different  measurement  geometry'  and  measurement  epoch,  GPS 
various  error  sources  which  include  ephemcris,  ionospheric  and  tropospheric  models,  receiver  coor¬ 
dinates  as  well  as  frequency  offset  between  the  reference  clocks  have  different  effect  on  GPS  single 
difference  observable. 
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Assuming  GPS  common  view  time  series  is  a  constant  linear  discrete  time  system.  For  a  specific  common 
view  time  series  i,  Z ,  ( k )  is  defined  as  the  optimal  estimate  of  the  i  “common  view  measuremen  t  X  ,  ( A: ) 
the  observation  equation  is: 


X((fc)  =  4>(^(fc)  +  f/ ,(*:) +  F,(fc)  (2) 

Where,  for  one  dimension  scalar  time  comparison,  the  coefficient  matrix  <f> ;  =  (  1 , 0 ) :  i  refers  to  the  i th 
common  view  measurement  (i  =  1  ,  rc  ) ,  k  refers  to  the  A:"'  day  (A:  =  J  is  the  systematic  part 

of  the  error  which  includes  the  frequency  offset  and  drift  between  the  reference  clocks,  propagation 
error,  ephemeris  error  as  well  as  receiver  coordinates  error.  £',-(/;) is  the  random  perturbation  part  of 
error  which  includes  the  random  perturbation  of  atomic  clocks  as  wel  I  as  the  measurement  noise  caused 
by  propagation,  ephemeris  and  receiver. 

The  systematic  bias  Ut( A;) between  the  i “time  series  and  the  j‘h  time  series  is  determined  by  the  GPS 
double  difference  observable  t  AB iy(t  (/)over  an  integral  data  length  (about  30  days)  and  the  average 
value  is  Bias  which  is  used  to  perform  correction  in  constant  relative  to  reference  time  scries.  The 
simulation  of  the  deterministic  properties  of  a  time  series  is  using  polynomial, 


X  t  =  a  +  bT  l  +  bT\  +  ...  +  bTnn  (3) 

Appjying  linear  regression  to  measurement  data  set  T k  and  X  t(k)  ,  (he  coefficients  of  the  polynomial 
are  determined  by  least  mean-square  estimate.  Thus  Z,-,  the  optimal  estimate  of  the  measurement  A'  „ 
is  obtained.  It  is  important  to  examine  the  confidence  in  the  result  by  taking  the  residual 
A(fc)  =  x f(k)  -  Z((k)  as  the  estimate  of  IF, (A:). 

First,  the  sampling  time,  x  ,  for  common  view  time  series  is  usually  one  day.  Assuming  F.  t  is  Gaussian 
noise  process,  it  includes  mainly  the  white  frequency  noise  and  dicker  frequency  noise.  So,  E ,  is  well 
modeled  by  equation  (4),  Allan  variance.  The  mean  squared  time  error  is  expressed  by  equation  (5). 

<i'(T)  =  (h0/2)x"‘  +  /<_,  21n(2)  (4) 

e2(x)  =  h0x  +  4/i_Jln(2)T2  (3) 


When  the  calibration  interval  xc  =  x  ==  1  day  ,  c2is  twice  as  large  as  u2,  and  it  is  not  divergent. 
Secondly,  the  typical  residuals  curve  for  a  common  view  data  length  of  about  135  days  is  in  Figure  1 


Let 


A  =  [(A(AO-Fm)/a]  =  300 


(6) 


F(.x)  =  V(2jO-V*Z/2  (7) 

where,  mean  value  F  m  =  0  so,  the  normalized  distribution  curve  (  F  ( x  )  -  X  ]  is  in  Figure  2.  It  conforms 
obviously  normal  distribution  N  (0 ,  a ). 

Hence  it  is  reasonable  to  have  a ,  the  standard  deviation  of  residuals  as  the  one  sigma  jitter  of  time 
series.  RMS,  the  standard  deviation  of  GPS  double  difference  gives  an  estimate  of  measurement  noise 
of  double  difference  observables.  The  variance  of  double  difference  between  the  i.  “  and  the  ;  common 
view  time  series  is  S  2.  Applying  the  N-corncr  hat  technique,  the  variance  estimate  RMS 2  of  the  i  "‘can 
be  completed: 

RMS2  =  (n-2)~'j^f  52-5^J  (8) 

6  =  (2n  -  4)" 1  (  J  J.92,](S2y  =  0)  (9) 

V  (*:-  1 )  y-  1  / 
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Thus,  RMS  can  be  considered  as  the  estimate  of  measurement  noise  V^. 

Bias,  Z  and  RMS  are  three  major  parameters  of  GPS  common  view  time  series.  The  common  view  data 
from  MJD4701 1  to  MJD47278  among  NRC,  NBS  and  USNO  are  divided  and  ealeu  la  ted  monthly.  Table 

1  gives  the  data  status.  Bias,  Z  and  RMS  are  calculated  in  each  data  section.  Taking  an  example,  Table 

2  lists  the  values  of  Bias,  Z  and  RMS  for  NRC-USNO  link  and  data  sections  from  1  to  5. 

The  results  in  Table  2  are  significant. 

1.  The  Biases  of  time  series  over  five  data  sections  are  quite  stable  the  variable  value  is  only 
several  nanoseconds.  By  using  linear  regression,  the  average  correlation  coefficients  are 
0.15  and  0.27  for  Z-Bias  and  RMS-Bias  respectively.  This  indicates  a  fact  that  systematic 
Bias  is  independent  of  random  perturbation,  hence  it  is  acceptable  to  correct  bias  in 
constant. 

2.  Usually  the  Z  value  is  larger  than  the  RMS  value.  The  £  value  seems  to  be  correlative  to 
the  data  section,  e.g.  the  Z  values  for  all  twelve  time  series  become  large  in  data  section 
3,  but  the  RMS  value  is  not  correlative  to  the  data  section.  These  conform  the  above 
analyses:  Z  is  the  sum  of  system  noise  and  measurement  noise,  while  RM.S  is  mainly  the 
measurement  noise. 

3.  The  values  of  Z  and  RMS  are  correlative  to  the  time  series  basically,  and  the  Bias  for  each 
time  series  is  quite  stable.  Thus,  the  lime  series  could  be  referred  to  as  independent  each 
other  and  to  weight  in  one  over  variance  is  optimum.  Because  the  correlation  coefficient 
of  Z-Bias  is  small  than  that  of  RMS-Bias,  to  have  Z  as  one  sigma  jitter  of  time  scries  is 
preferred. 

4.  The  above  results  are  obvious  not  only  for  N RC-USNO  link  but  also  for  NRC-NBS  link. 

UNBIASED  WEIGHTING  METHOD 

Based  on  the  above  analysis,  an  algorithm  for  processing  data  of  GPS  common  view  time  comparison 
is  recommended  to  be  referred  to  as  Unbiased  Weighting  Method  fUBWA).  It  has  the  following  fea¬ 
tures. 

1.  Averaging  the  double  difference  observables  during  an  integral  length  30  days’,,  the 
systematic  biases  existing  between  every  two  time  scries  are  corrected  in  constant  monthly 
relative  to  a  reference  time  series. 

2.  The  one  sigma  jitter  of  each  time  series  is  determined  by  the  standard  deviation  Z  of 
residuals  of  polynomial  regression,  which  defines  weights  in  one  over  Z  square. 

3.  Using  raw  measurement  data  at  all.  no  interpolation  or  extrapolation  for  losing  data  points. 
For  bad  points,  it  can  be  identified  and  removed  by  computer  program. 

The  steps  and  formulae  for  UBWA  are  as  follows.  The  ktfl  day  and  the  ilh  path  common  view  data 
between  two  stations  is  Xj^j,  and  the  i^1  path  common  view  data  form  a  time  series: 

X ;  ( A: )  (fc-l.-'H)  (  i  0  ) 

Where,  m  is  total  measurement  days,  in  practice  M=30  as  a  data  section,  i-..  n,  n  is  the  number  of 
common  view  paths  each  day  between  two  stations.  So  total  has  n  ume  series.  Do  least  square  t  it  for 
the  i*h  time  series  and  calculate  the  standard  deviation  Z ,  of  residuals,  then  the  optimum  weight  for  the 
ith  series  is: 

t/,=  l/Zf  (ii) 

Choose  a  reliable  and  continuous  time  series  as  reference  one,  namely  the  jta.  The  double  differences 
for  all  other  series  relative  to  reference  onexan  be  calculated, 
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(12) 


/' tj(k')  =  X f(fc)  -  X j(k)  (i  =  1  ,  n  -  1  ;  j  ¥=  i;  fc  =  1  ,  m’) 

The  average  systematic  bias  Sj  between  the  qh  series  and  series  is: 

m 

St  =  Y.  Ti/W/rn 

k-  1 

The  Sj  value  is  used  to  correct  the  daily  data  of  the  i*h  time  scries, 

Xt^X,W*St  (j*  j;k<=  1  ,m) 

Xik  =  Xt(k)  ( t  =  j)k=  1  ,  m ) 


(13) 


(14) 


Considering  the  losing  data  points,  the  daily  weight  for  the  i^1  time  series  needs  to  be  calculated, 


h/«,  =  l/i 

A7(fc)7  0 

IA,  =  0 

o 

II 

} 

(1.5) 

So, 

n 

AWtk  =  Wlk/Y.Wik 

(16) 

t“  1 


Finally,  the  weighted  average  of  n  time  scries  at  the  measurement  epoch  of  the  j^1  time  series  can  be 
obtained, 

n 

x k=  Y,  AW.ikXik  (A:  =  1  ,  in)  (17) 

1 


RESULTS 

The  common  view  data  among  NRC,  NBS  and  USNO  are  computed  by  using  three  different  methods 
from  data  files  extending  for  a  length  of  about  nine  months.  The  time  series  of  final  weighting  average 
arecompared  in  three  performance  characteristics.  The  three  met  hods  are:  Unbiased  Weighting  method 
(UBWA),  Direct  Weighting  method  (DWA)  which  does  not  correct  the  systemat  ic  bias,  as  well  as  Equal 
Weights  method  (EWA).  The  three  characteristics  are:  the  standard  deviation  I  of  residuals  of  poly¬ 
nomial  regression,  the  frequency  stability  Allan  variances  at  different  sampling  time  and  the  time  dis¬ 
continuity  caused  by  increasing  or  decreasing  the  number  of  time  series.  Results  showed  that  UBWA 
is  prior  to  the  others. 

1.  Data  file  is  divided  into  ten  time  sections  showed  in  Table  1  and  computed  in  three 
methods  monthly.  The  weighted  average  scries  UBWA,  DWA  and  EWA  are 
obtained  respectively.  The  standard  deviations  2  of  residuals  of  three  time  scries 
are  listed  in  Table  3.  The  UBWA  results  are  better  than  the  others,  they  are  about 
4.9  ns  for  NRC-NBS  link  and  about  2.4  ns  for  NRC-USNO  link.  The  lvalue  is  the 
indication  of  time  perturbation  of  time  series. 

2.  The  frequency  stability  Allan  variances  at  sampling  times  of  1,  2,  3,  7,  12,  16,  and 
20  sidereal  days  are  calculated  for  UBWA,  DWA  and  EWA  weighted  average  lime 
scries  respectively.  Results  are  plotted  in  Figure  .3  which  has  coordinates  log 
logUy(t)-  10 log's-  and  x  in  seconds.  UBWA  scries  has  the  best  frequency 
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stabilities.  The  a  ( x  )  is  better  than  1  part  in  10^  for  sampling  times  of  five  days 
and  longer,  this  meets  the  requirements  of  high  precision  time  and  frequency 
comparison  between  timing  laboratories. 

3.  Because  of  the  existence  of  systematic  bias,  the  time  discontinuity  could  appear  at 
the  conjunction  point  between  two  consecutive  data  section  when  the  weights  of 
time  series  changed  largely  or  the  number  of  time  series  increased  or  decreased 
especially.  Due  to  correction  of  Bias,  this  effect  should  be  not.  obvious  for  UBWA 
method.  During  data  sections  1  to  5  from  MJD  47144  to  47278,  the  number  of  time 
series  was  changed  intentionally  to  exam  the  possible  time  discontinuity.  The 
computed  results  are  listed  in  Table  4  and  the  time  discontinuity  is  the  smallest  for 
UBWA  method. 


4,  In  UBWA  method,  it  is  required  to  choose  a  reference  time  series  and  to  correct 
systematic  biases  for  other  time  series  relative  to  the  reference  one.  For  large 
geographic  area  or  multistation  GPS  common  view  data,  it  is  not  possible  Lo  find  a 
common  reference  time  series.  But,  this  can  be  solved  by  using  reference  time  series 
relay  method,  and  time  accuracy  could  be  ensured. 
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Table  1 .  Data  status 


section 

MJD 

Interval 

(days') 

No.  of  time 
scries 

1 

47248-4727S 

31 

12 

2 

47218-47248 

31 

12 

n 

J) 

47188-47218 

31 

12 

4 

47158-47188 

31 

12 

5 

47144-47158 

15 

12 

6 

47113-47143 

31 

10 

7 

47083-47113 

31 

10 

8 

47053-47083 

*->  -i 

1 

10 

9 

47023-47053 

31 

10 

10 

47011-47023 

13 

10 
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Table  2.  Parameter  estimate  for  common  view  time  series  (NRC-USNO)  (ns) 


Path 

031 

061 

091 

092 

093 

111 

112 

113 

121 

122 

123 

131 

1 

-4 

23 

-13 

14 

-1 

17 

23 

20 

-18 

20 

16 

2 

-5 

24 

-13 

-12 

-4 

-16 

22 

18 

-16 

20 

20 

Bias 

3 

-1 

28 

-7 

-13 

-4 

-15 

22 

24 

-18 

25 

16 

4 

-2 

24 

-13 

-14 

-5 

-16 

23 

21 

-19 

22 

17 

5 

3 

22 

-13 

-7 

0 

-14 

25 

21 

-16 

21 

19 

1 

3.6 

3.6 

3.4 

4.1 

4.6 

3.4 

3.8 

3.4 

4.6 

4.5 

5.1 

3.7 

2 

3.8 

2.5 

3.6 

3.1 

2.6 

4.6 

3.3 

3.4 

4.2 

3.3 

4.9 

3.3 

>: 

3 

8.6 

8.4 

12.5 

7.4 

8.9 

8.0 

9.9 

8.6 

9.4 

9.2 

10.2 

6.9 

4 

4.8 

5.4 

5.2 

5.2 

4.9 

4.5 

5.3 

4.0 

4.8 

3.4 

5.3 

5.6 

5 

5.2 

4.0 

4.9 

4.2 

5.8 

4.5 

4.7 

6.0 

6.1 

3.3 

6.1 

4.8 

1 

2.8 

3.6 

3.3 

3.2 

4.2 

4.0 

3.0 

1.6 

4.4 

3.5 

4.2 

3.4 

2 

3.1 

2.1 

5.2 

4.5 

2.8 

4.4 

3.0 

2.7 

3.6 

2.7 

4.7 

3.0 

RMS 

3 

2.7 

2.1 

5.8 

2.3 

2.6 

2.2 

4.8 

2.9 

2.0 

3.3 

5.0 

3.8 

4 

3.2 

3.2 

4.4 

2.8 

2.5 

1.9 

4.6 

2.6 

5.6 

2.4 

4.5 

3.0 

5 

1.6 

2.6 

5.8 

4.9 

7.8 

1.5 

1.7 

2.1 

4.8 

2.2 

3.0 

3.4 

Table  3. 1  values  of  weighted  average  time  scries 


Data  section 

5  ->  4 

4  ->  3 

3  ->  2 

2  ->  1 

Number  of  time 

12 ->  10 

10->  12 

12 ->  10 

10 ->  12 

series 

Conjunction 

47158 

47188 

47218 

47248 

Date 

NRC  UBWA 

+  1.7 

-2.0 

+  1.9 

-2.0 

/  DWA 

-1.7 

-7.3 

+2.3 

+3.2 

NBS  EWA 

+4.2 

-4.2 

+4.0 

-3.7 

NRC  UBWA 

-2.0 

+2.4 

-0.4 

-1.1 

/  DWA 

2.2 

-1.2 

+5.4 

-5.2 

USNO  EWA 

+3.6 

-2.3 

+3.5 

0.7 

524 


RESIDUALS 


Figure  2.  Normal  distribution  of  residuals  Figure  3.  Allan  variance  results  for  UBWA,  DWA  and  EW'A 


